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We examine the possibility that thermal, rather than quantum, fluctuations are responsible for 
seeding the structure of our universe. We find that while the thermalization condition leads to 
nearly Gaussian statistics, a Harrisson-Zeldovich spectrum for the primordial fluctuations can only 
be achieved in very special circumstances. These depend on whether the universe gets hotter or 
colder in time, while the modes are leaving the horizon. In the latter case we find a no-go theorem 
which can only be avoided if the fundamental degrees of freedom are not particle-like, such as in 
string gases near the Hagedorn phase transition. The former case is less forbidding, and we suggest 
two potentially successful "warming universe" scenarios. One makes use of the Phoenix universe, 
the other of "phantom" matter. 



I. INTRODUCTION 

A major success of the Big Bang theory is its abil- 
ity to account for the detailed structure of the Universe, 
such as galaxy clustering, or the temperature fluctua- 
tions in the cosmic microwave background. Paramount 
to this picture is the phenomenon of gravitational in- 
stability, by means of which primordial small departures 
from homogeneity can grow into the observed structures. 
Within the Big Bang theory the required primordial fluc- 
tuations are treated merely as initial conditions. Vacuum 
quantum fluctuations in inflationary scenarios have been 
shown to lead to the required initial conditions, fitting 
current large scale structure data. It is, however, impor- 
tant to know just how unavoidable is the conclusion that 
cosmic structures have a quantum origin. Could these 
primordial fluctuations have a thermal origin instead? 

The possibility that we are descendent from thermal 
fluctuations was advanced by Peebles, in his book Prin- 
ciples of Physical Cosmology, pp. 371-373 [1]. Peebles 
pointed out that if the Universe was in thermal equilib- 
rium on the comoving scale of 10 Mpc when its temper- 
ature was T = 10 11 Gev, then the observed value of uio 
could be explained. This fascinating remark leaves sev- 
eral questions unanswered. Such a scenario may explain 
the observed value of oio, but what about fluctuations 
on other scales? Furthermore, thermal fluctuations are 
not strictly Gaussian - does this scenario conflict with 
observations? 

Motivated by these unsolved issues, in this paper we 
go further and examine under which conditions primor- 
dial thermal fluctuations lead to a Harrison- Zeldovich 
(HZ) [2] spectrum of approximately Gaussian fluctua- 
tions. Such a scenario would fit all existing data in the 
same way that the usual inflationary quantum fluctua- 
tion scenario does. Indeed, the only potentially distin- 
guishing feature would be a different signature (or ab- 
sence thereof) of gravitational waves (tensor modes) in 
the thermal scenario. 

A matter not addressed by Peebles is how to es- 
tablish thermalization on the relevant scales. In pure 



Big Bang cosmology all observed fluctuations spanned 
causally disconnected regions when the universe was at 
the required high temperatures to induce the appropri- 
ate level of fluctuations. In addition, simply postulating 
exact thermalization over all scales, say at Planck time, 
leads to grossly inappropriate results (thermal fluctua- 
tions are white-noise rather than scale invariant). Hence, 
we should assume that exact thermalization only applies 
to sub-horizon modes, and that there is a mechanism for 
pushing sub-horizon thermal modes outside the horizon, 
where they freeze and become non-thermal. Since differ- 
ent modes freeze at different temperatures, it may be that 
the final super-horizon spectrum is indeed of HZ type. 

We shall consider three types of mechanisms for push- 
ing modes outside the horizon: accelerated expansion, 
varying speed of light (VSL) and a contracting universe. 
In the first case we consider inflation models driven by 
a dominant component of thermal radiation, as in [3], 
where deformed dispersion relations affect common ra- 
diation at high temperatures. This is not to be con- 
fused with inflationary models in which there is a finite 
radiation component during inflation [4,5]. There the 
dominant component is always the inflaton field (even 
though the thermal bath drives inflaton fluctuations). 
The second possibility is a varying speed of light (VSL), 
either in the form of a space-time field (c(x^)) [6-9], 
or as an energy-dependent effect (c(E)) [11]. There 
are VSL models [8] in which quantum vacuum fluctu- 
ations can produce the HZ spectrum [10]. However, in a 
large class of VSL scenarios the universe is never vacuum 
dominated. Hence subhorizon scale thermalization is the 
reason for the apparent homogeneity of the Universe - 
and likewise thermal fluctuations are responsible for the 
primordial fluctuations. In the third case, we consider 
the possibility of thermal fluctuations seeding structures 
within Lemaitre's Phoenix universe [12]. 

The paper is organized as follows. In Section II we ex- 
amine the statistics of thermal fluctuations. In Section 

III we derive the necessary conditions for thermal fluc- 
tuations to have a scale-invariant spectrum. In Section 

IV we show that these conditions cannot be fulfilled in 
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universes dominated by radiation comprised of conven- 
tional particles and in which the temperature decreases 
with time. At the end of Section IV and in Section V 
we propose a set of models that may bypass this no-go 
theorem. We summarize our results in Section VI. 



II. THE STATISTICS OF THERMAL 
FLUCTUATIONS 

Most of the currently available data is consistent with 
the hypothesis that primordial fluctuations are Gaussian 
distributed [13]. A possible exception is the report of 
non- vanishing "inter-scale" components of the CMB bis- 
pectrum calculated from the four year COBE-DMR data 
[14]. New CMB measurements, in particular by the MAP 
satellite, are expected to provide tighter constraints on 
the amount of cosmological non-Gaussianity. 

Based on current observational evidence, one has to 
require that at least on very large scales, primordial fluc- 
tuations must be sufficiently well described by a Gaussian 
distribution. Here we encounter the first obstacle, since 
strictly speaking thermal fluctuations are not Gaussian. 
In what follows we shall show that these fluctuations are 
Gaussian to a very good approximation under the same 
set of conditions which assure thermalization. 

Fluctuations in a thermal (canonical) ensemble can be 
determined from the partition function 
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where (3 = T . This is true even if deformed dispersion 
relations are introduced [11]. Hence the total energy U 
inside a volume V is given by: 



U=(E)= E P e-^ —W 
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In general this integral needs not be proportional to T 4 , 
and indeed under deformed dispersion relations U cx T 7 , 
with 1 < 7 < 4. However, because energy is an exten- 
sive quantity, we always have that U = pV, that is U is 
proportional to the volume. If 7 ^ 4 there is a preferred 
length scale It, and we can choose units so that 
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(specifically we use units such that G = U = Co = ks = 1, 
where c is current value of c, and neglect factors of order 
!)• 

The energy variance is given by 
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and so the relative variance is 
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where prime denotes differentiation with respect to tem- 
perature. We see that whereas the amplitude of these 
fluctuations is model dependent, their spectrum (i.e. the 
fact that cry cx 1/V) is not. The white noise nature of 
thermal fluctuations follows from the fact that energy is 
an extensive quantity (i.e. proportional to the volume). 

The fact that afj cx 1/V leads to an interesting heuris- 
tic interpretation of thermal fluctuations. It seems to im- 
ply that thermal fluctuations may be seen as a Poisson 
process involving a set of regions with coherence length A 
dependent only on the temperature (and not the sample 
volume). In a volume V there are n = V/X 3 such regions, 
so that a Poisson process results in variance 
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implying that any white noise spectrum of fluctuations 
may be seen as a Poisson process. The dependence A(T) 
can be inferred from eq. (5) and in general takes the form 
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translating, for (3), into 
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We see that A cx T^ 1 only for 7 = 4. For 7 = 1 (re- 
alized in non-commutative geometry [11]) A is temper- 
ature independent and equals the length scale of non- 
commutativity. In general the thermal coherence length 
decreases with increasing temperature, except if 7 < 1. In 
the exceptional case 7 < 1 the coherence length increases 
with the temperature and the relative energy fluctuations 
anomalously increase with the temperature; in Section IV 
we shall rule out this exceptional case. 

The non-Gaussianity of thermal fluctuations may now 
be studied in terms of the cumulants of the distribution. 
The third order centered cumulant is given by [15] 
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so that the relative skewness is 
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Hence for large volumes S3 C 1. Likewise for higher 
cumulants 
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and so 
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with the proportionality constant c„(7) = 7(7 + 1)...(7 + 
n — 2)/7™/ 2 . Setting V = L 3 we find the unsurprising 
result that this can be written as 



Cn(7) ( x 



7 5 



3(f-l) 



(13) 



and so if the volume under study is much larger than 
the thermal coherence length as defined above we have 
indeed that s n -C 1. 

We conclude that thermal fluctuations are Gaussian to 
a very good approximation if, when they leave the hori- 
zon, the thermal coherence length is much smaller than 
the horizon size. The departures from Gaussianity always 
lead to positive cumulants and these decay with the order 
n as in eq. (13). Since the condition for thermalization is 
that the coherence length is much smaller than the scales 
under study, we may conclude that Gaussianity is part 
and parcel of the self-consistency conditions for studying 
thermal fluctuations in thermal equilibrium. 



outside the horizon and are allowed to re-enter it at a 
later time 1 . Then, at first horizon crossing we have 
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where k are comoving wave numbers, a is the scale factor 
and c is the speed of light 2 . The modulus sign in eq. (17) 
accounts for the possibility of such a crossing happening 
during a contracting phase in a bouncing universe - a 
model discussed in subsection VB. We parameterize the 
temperature dependence of the first horizon crossing with 

k h = T^r 1 , (is) 

where Ai is some characteristic length scale. We will 
consider two types of scenarios: those in which the tem- 
perature of the radiation is decreasing with the evolution, 
and those in which the temperature is increasing. 

If the universe cools as it expands then solving the hori- 
zon problem and the existence of a first crossing require 
that p < 0. More generally, one needs dkh/dT < or, 
from eq. (18), 



III. THE POWER SPECTRUM 



Gaussian fluctuations are fully described by their two- 
point function. This can be encoded in afj, but is more 
often expressed in terms of the power spectrum P(k) = 
(\Sk\ 2 ), where Sk are the Fourier modes of the density 
contrast Sp/p. The two can be related via the integral 
[16]: 
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where Wp(kL) is a filter function and L <~ V 1 ^ 3 is the 
smoothing scale (so that U — pV). Assuming a power 
law dependence for the power spectrum, P{k) = A 2 k n , 
and, for instance, a Gaussian filter function, Wp(kL) = 
e -k l /2^ integration of eq. (14) gives 
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where A 2 = A 2 Y (2±3) /(4tt 2 ), that is, IrjA. By com- 
paring cqns. (5) and (15) one can see that thermal den- 
sity fluctuations have a white noise spectrum (n = 0) 
with amplitude: 
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where we have ignored factors of order 1. This result 
only applies to modes which are in causal contact, i.e. 
sub-horizon modes. More precisely it only applies when 
self-gravity is negligible, and so to modes smaller than 
the Jeans length. 

Suppose we have a model in which a certain range of 
Fourier modes of thermal density fluctuations are forced 
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Hence, while modes are being forced outside the horizon 
one must have either accelerated expansion or a decreas- 
ing speed of light, or a combination of both. 

If the universe heats up as it evolves in time, then p > 
(or, dkh/dT > 0). We consider this case in some detail 
in Section V. 

In both types of scenarios we will be interested in iden- 
tifying conditions for a HZ spectrum of density fluctua- 
tions to be left outside the horizon. If there is no sig- 
nificant evolution of the gravitational potential <fi outside 
the horizon while the modes are being pushed out (this is 
usually enforced by requiring that the equation of state 
remains more or less constant), then density fluctuations 
have a HZ spectrum when the equal-time power spectrum 
of </> has a form 

k 3 P^ = B 2 , (20) 

with B :=y 10~ 5 . If in addition at horizon crossing we 
have 4>h w Sh, then we need 3 : 



1 This condition is more restrictive than requiring that the 
horizon problem be solved - for instance the Milne universe 
(with a oc t) does not have horizons, and yet sub-horizon 
modes are not pushed outside the horizon. The reader is also 
referred to Ref. [17] for a simple condition for solving the 
horizon and flatness problems. 

2 c is the speed of light at the time of the horizon crossing, 
which could be different from Co 

a This may be seen as an independent assumption, or justi- 
fied using Einstein's equations. The perturbed Friedmann 
equations in the comoving longitudinal gauge imply that 
k 2 (f> ~ a 2 p8 and, using the Friedmann equation, a 2 
and eq. (17), we obtain 4> « (kh/k) 2 8. Hence 4>h ~ Sh- 
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which is sometimes used as the definition of the HZ spec- 
trum. Since fa ~ 5u , this expression is true in any gauge 
because the Sh defined in the various gauges are all pro- 
portional to each other. 

If we identify these fluctuations with thermal fluctu- 
ations about to leave the horizon, then, using eq. (16), 
this is equivalent to 
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(where again we have neglected factors of order 1). eq. (3) 
then leads to 
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Taking the time derivative of the above expression leads 
to 



(where w is the equation of state) , and if we parameterize 
the dependence of the speed of light on the temperature 
as: 



(29) 



Then from eqns. (17), (18), (28) and (29) the condition 
for scale invariance reduces to: 



7 = 2 1 + w 
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IV. A NO-GO THEOREM FOR COOLING 
UNIVERSES 
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Since kh > is the condition for modes to be leaving the 
horizon, it follows from eq. (24) that 
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Hence scale invariance requires that 7 < 1 or 7 > 1 
depending on whether the Universe is getting colder or 
hotter. We recall that 7 < 1 is equivalent to saying 
that the thermal coherence length A increases with the 
temperature (see eq. (8)), or that the fractional energy 
fluctuations increase with the temperature. This is quite 
anomalous and we shall rule it out explicitly in the next 
section. Hence we are left with warming universes as a 
possibility for HZ fluctuations of thermal origin. 

These results may be expressed more quantitatively by 
noting that Eqs. (18) and (23) lead to the condition for 
scale invariance: 



7 = 1 + 3m , 
or more generally the expression for the tilt 
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In our argument so far we have abstained from using 
Einstein's equations. Instead, we have treated eq (21) 
as an independent assumption and also assumed that the 
gravitational potential would freeze outside the horizon. 
However we can go further if we are prepared to use Fried- 
mann equations: 
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If at all times the universe has been cooling (T < 0), 
then we arrive at a forbidding set of conditions for a scale- 
invariant spectrum. In this case kh must decrease with 
the temperature (fi < 0). As already noted, from (25) it 
then follows that this requires 



7 < 1 
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This condition is a strict inequality: 7 = 1 leads to either 
n = 4 (if the modes are indeed being pushed out of the 
horizon), or to n = (if the Hubble length stagnates). 

If 7 < 1 the fractional amplitude of thermal fluctu- 
ations anomalously increases with the temperature. As 
equation (8) shows, this also implies that the thermal cor- 
relation length increases with the temperature. Below we 
present a no-go theorem which shows that it is unlikely 
that this condition is satisfied assuming that the weak 
energy condition is satisfied and that the fundamental 
degrees of freedom are particle-like. 

Consider radiation in thermal equilibrium with a cer- 
tain dispersion relation p(E) which becomes the usual 
p 2 = E 2 at sufficiently low energies. The energy density 
is proportional to the integral: 



p(T) oc J(T) = J dE 
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where the integration is over all allowed values of energy. 
For convenience, let us define F(E) = Ep 2 (E) \dp/ dE\ 
and re-write J(T) as: 



(34) 



In order to have 7 < 1, f(T) must be a decreasing func- 
tion of temperature at sufficiently high values of T, i. c. 
f'(T) < . Let us evaluate f'(T) : 
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If we only consider non-negative energies, then the factor 
under the integral appearing in front of the square brack- 
ets is non-negative, while the expression inside the brack- 
ets is a non-negative, monotonically increasing function 
oiE/T. Hence, f'(T) > for all T. Thus, the inequality 
(32) cannot be satisfied. 

The above proof is quite general and is valid for all 
models that would aim to achieve 7 < 1 by modify- 
ing the dispersion relations without altering the statis- 
tical properties of the gas. In particular, this proof im- 
plies that the modified dispersion relations considered in 
Refs. [18,11,3]) could not result in a HZ spectrum. 

The no-go argument assumed that the thermalized 
radiation was made of particles obeying Bose-Einstcin 
statistics. A radically different statistics is likely to be 
needed in order to obtain the desired exponent in the 
Stcfan-Boltzmann relation. Thus the only way we see of 
bypassing this no-go theorem is to allow for non-particle 
like degrees of freedom. We conclude this section by sug- 
gesting a scenario which may make use of this possibility. 



A. Saturating temperature 

An example of a system in which it is possible to have 
7 < 1 is a gas of strings at temperatures close to the so- 
called Hagedorn temperature, Th [19]. In a gas of strings, 
the number of degenerate states increases exponentially 
with energy [20] and the canonical partition function di- 
verges for all T > Th- This does not necessarily mean 
that temperature higher than Th are unphysical. In fact, 
all physical quantities, such as energy density and spe- 
cific heat are actually finite at T > Th [21]. In [22] it 
was suggested that Th corresponds to a phase transi- 
tion, somewhat analogous to the deconfining transition 
in QCD. At temperatures close to Th, the canonical en- 
semble description of string gases becomes invalid due 
to increasingly large energy fluctuations [23]. One must 
use the microcanonical ensemble instead, which is well- 
defined only if all spatial dimensions were compactified 
[21] 4 - 

At least within the canonical ensemble formalism, Th 
can be interpreted as the limiting temperature of the gas 
- as energy is increased, the temperature remains con- 
stant. In the language of eq. (3) this corresponds to 
7 — > 0, in agreement with the constraint (32). A straight- 
forward examination of eq. (30) with 7 = shows that 
scale- invariance can be satisfied if cither w = — 1, as in 
inflation, or if a = 1/3, as in VSL models. 



4 In Ref. [21] it is further suggested that in this picture one 
needs a mechanism which would later make three of the spa- 
tial dimensions sufficiently large for us to live in. 



String-driven inflationary models, making use of the 
existence of a limiting temperature, have been considered 
in the late 1980's [24,25]. More recently, in Ref. [26], it 
was proposed that winding modes of open strings on D- 
branes above the Hagedorn phase transition can provide 
the negative pressure necessary to drive inflation. In par- 
ticular, it was suggested that one could achieve a period 
of exponential inflation (with w—-l) if not all transverse 
dimensions supported winding modes [26]. 

It would be very interesting to investigate if thermal 
fluctuations could indeed be viable candidates for struc- 
ture formation in these models and the degree of fine- 
tuning it would involve. 

The other possibility, that of a VSL theory with a = 
1/3, is currently lacking a specific model realization. 



V. WARMING UNIVERSES 

But it could be that at the early stage when the modes 
are being pushed out of the horizon the universe is get- 
ting hotter. Such is the case of thermal radiation with 
7(1 + w) < 0. If 7 > 0, denser radiation means hotter 
radiation; however if w < —1 the universe gets denser 
(and so hotter) as it expands. Alternatively we could 
have w > —1, so that the universe gets less dense as it 
expands; but then if 7 < this translates into a higher 
temperature. 

Another possibility is a stage of radiation injection, 
either from particle/antiparticle annihilation, from false 
vacuum decay, or from a cosmological constant discharge. 

Yet another possibility is the Phoenix universe of 
Lemaitre [12], where modes would be pushed outside the 
horizon with temperature increasing during the contract- 
ing phase. 

If the universe gets hotter in time, we need kh to in- 
crease with time, and with temperature, so that fi > 0. 
A necessary condition for scale invariance is then 



7 > 1 



(36) 



bypassig the no-go theorem in the previous section. 
Again, from eqns. (17), (18), (28) and (29) we obtain 



= 7 

or equivalently 



1 



1 



2 3(1 +w) 



a 1 + 3w 
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a>0 



(37) 



(38) 



Conditions (30) and (31) still apply. 

We now consider particular solutions to these condi- 
tions. 
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A. A phantom phase 

"Phantom" matter [27] exhibits an equation of state 
with w < — 1, and it may constitute the dark energy of 
the universe. It has also been conjectured that normal ra- 
diation at high temperatures could behave like phantom 
matter [3]. For these models there is a critical density, 
p c , such that w > — 1 for p < p c , while for p > p c one 
has w < — 1. If the Universe starts off with p > p c and 

2 

expanding then a oc (— t) 3 ( 1 +™> . As the universe expands 
it gets denser and hotter. Eventually a phase transition 
brings it to the sub-Planckian regime. 

In such a scenario there is hyper-inflation, so the modes 
are pushed out of the horizon without a VSL. However, 
in order for density fluctuations to have a scale invariant 
spectrum, the condition in eq. (30) must be satisfied. 
Since 7 > 1 and w < —1, we find that 



1 

a> - 



1 



1 

> - . 
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Hence, in order to obtain scale-invariance, one needs a 
VSL. 

Regarding the spectrum's amplitude, from eq. (23) one 
can obtain the requirement 



X 1 » w B 2 w 1CT 10 



or, using eq. (26) to express p in terms of A, 
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Thus, given a value of 7, eq. (41) constrains the ratio of 
the two characteristic length scales It and Ai. Within 
any specific model providing a relation between the scale 
factor a and the temperature T it should be possible to 
express Ai in terms of It- Consequently, these two length 
scales should not be considered as independent. 

We may expect one of the models based upon the for- 
malism presented in [11] to satisfy these conditions, plac- 
ing further restrictions upon deformations of dispersion 
relations. The spectrum's amplitude in this case will be 
fixed by the ratio between the non-commutative length 
scale (related to It) and the Planck length (related to 
Ai). 



B. A bouncing universe 

In the bouncing model a closed universe goes through 
a series of cycles starting with a Big Bang and expansion, 
followed by re-contraction and a Big Crunch. Every time 
the universe approaches the crunch it bounces into a new 
bang, and if entropy increases at the bounce the new 
cycle lasts longer (i.e. expands to a larger and colder 
turnaround point) in each new cycle. Although the idea 



goes back to Lemaitre [12], it has often been forgotten 
and revived a few times, e.g. more recently in [28]. 

It is questionable that this model dispenses with fine- 
tuning, specially when dealing with the homogeneity and 
entropy problems. It may be that just too much junk (en- 
tropy) is left over from previous cycles. This can usually 
be avoided combining the bouncing universe with infla- 
tion, such as in its ekpyrotic incarnation [29]. 

Another possibility is that the speed of light decreases 
at the bounce, so that all the relevant scales we see to- 
day were sub-horizon modes at some point deep in the 
radiation epoch in the contracting phase. Hence homo- 
geneity (and absence of pervading coalescing black holes) 
could have been established on all the relevant scales 
for the current cycle. But more importantly in this sce- 
nario thermal fluctuations could become the Harrisson- 
Zeldovich spectrum of initial conditions which we ob- 
served in our current cycle. We now illustrate how this 
could happen. 

As long as 10 > —1 (and 7 > 0) the universe heats 
up in time during the contracting phase. Hence a nec- 
essary condition for scale invariance is that 7 > 1, evad- 
ing the no-go theorem presented in Section IV. Modes 
must leave the horizon as the universe contracts, and as 
eq. (37) shows this can be achieved without a VSL if 
w > —1/3, i.e. the universe undergoes the necessary ac- 
celerated contraction when the strong energy condition 
is satisfied. A general condition for scale invariance is 
then eq. (30). We find the notable result that standard 
radiation (w = 1/3, 7 = 4, a = 0) satisfies the condi- 
tion for scale invariance. Hence, as long as the scales we 
can observe today were sub-horizon deep in the radia- 
tion epoch during the contracting phase of the previous 
cycle, thermal fluctuations in standard radiation lead to 
a Harrisson-Zeldovich spectrum of initial fluctuations in 
the new phase 5 . 

A problem arises with the amplitude, which may be 
remedied if there is a drop in the speed of light at the 
bounce. For 7 = 4 the scale It drops out of the problem 
(cf. eq. (3)). Also for w = 1/3 and a = we have 
kh oc T, so p = 1 and the scale Ai disappears from 
eq. (18). Hence the ratio of these parameters disappears 
from the expression for the amplitude, as can be seen 
directly from Eqn. (41); the amplitude becomes then a 
combination of numerical factors of order 1 clearly in con- 
tradiction with observations. This can also be illustrated 
using the 10 Mpc comoving scale as the normalization 
point, as in Peebles' example [1]. Rephrasing Peebles 
argument in terms of this model, we know that in the 
case of a perfectly symmetric bounce the right normal- 
ization 010 can be obtained if the 10 Mpc comoving scale 



5 This assumes certain mode matching at the bounce. In 
general, the picture will be much more complicated. This 
issue is discussed in subsection V B 1 
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between time and temperature is symmetrical around the 
bounce, but for c_/c + <C 1 the 10 Mpc comoving scale 
would leave the horizon in the contracting phase much 
later and at a higher temperature, as illustrated in Fig. 1. 

Using t/(lsec) w (T/(lMcv))- 2 , we find that U = 
1CT 28 sec before the bounce the temperature is 10 11 Gev. 
The comoving horizon at this time is 10 Mpc across if 



C-U = lOMpc— 

a 



(42) 



10 



21 



(43) 



Since a*/a = Tc7Ms/10 1:L Gev (where T C mb ~ 2.7 K) 
we get 

c_ 
c+ 

Note that unlike other VSL arguments, which lead to 
lower bounds on c_/c + , this argument leads to an iden- 
tity: the spectrum amplitude results directly from a given 
value of c_ /c+ . It may be possible to obtain a similar ef- 
fect by another way of making the bounce asymmetrical. 
For instance if G and h vary differently at the bounce a 
different constraint is obtained. 

Certainly, the validity of the above discussion depends 
upon the matching of growing and decaying modes dur- 
ing the bounce dynamics - something which is unknown 
in the absence of a specific model for the bounce. We 
conclude with a discussion of the two main uncertainties 
faced by this model. 



FIG. 1. The normalization problem in bouncing universes 
and a way to solve it. In the case of a perfectly symmetrical 
bounce, density perturbations on 10 Mpc scale exit the hori- 
zon during the contracting phase at temperatures much lower 
than 10 11 GeV. The problem could be fixed if the speed of 
light in the contracting phase (c_ ) had a much bigger constant 
value than it has today (i. e. c_ S> c+). 

leaves the horizon in the contracting phase when the uni- 
verse was at 10 11 Gev. However, if the relation between 
temperature and time on cither side of the bounce is 
symmetric, this comoving scale would leave the horizon 
much earlier in the contracting phase, when the universe 
was much colder and therefore the fractional fluctuations 
were much larger. 

The problem with the normalization mentioned above, 
could be resolved if the bounce was asymmetrical, allow- 
ing for the 10 Mpc comoving scale to leave the horizon 
when the universe was indeed at 10 11 Gev. One possi- 
ble scenario is when the speed of light is constant during 
the previous (c = c_) and current (c = c+) cycles, but 
c_ 3> c + , i.e. the speed of light decreases at the bounce. 
One can make additional simplifying assumptions just to 
illustrate the point. Let us assume that the fundamental 
constants G and h also change at the bounce in such a 
way so as to preserve the relation between the temper- 
ature and time (this can be achieved if h 3 c 5 /G is the 
same before and after the bounce). Hence the relation 



1. Mode matching at the bounce 

Bounce dynamics will in general mix the modes in ci- 
ther phase [34,35], a matter which may upset our previ- 
ous arguments. In the contracting phase we have 

= A^ + ^ji (44) 

and in the expanding phase 

= A + cj> + ^ (45) 

where <fr is the gravitational potential and t = at the 
bounce. Matching these modes depends on the bounce 
dynamics, and everything we have said previously as- 
sumed that _B_ mode was not excited and that the A_ 
mode was matched onto the A + mode. 

If the B- mode is excited, but then matches on to the 
B + mode, then the earlier discussion still applies. How- 
ever if the £?_ mode matches onto the A + mode then the 
conclusions in this paper have to be revised: in addition 
to the HZ spectrum predicted here, there would be a very 
red component (n = —3). Another pathological case is 
when A_ and B + are excited, while £?_ = A + = 0. 

We defer to a future publication a more careful study 
of this situation. It certainly depends on how the bounce 
is actually produced. Mode matching in pre-big-bang 
[36] and ekpyrotic [29] cosmologies has been discussed in 
considerable detail in Rcfs. [34,35] 
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2. Entropy production at the bounce 

Another matter which may change the spectrum and 
normalization of the fluctuations is entropy production 
at the bounce (which almost certainly occurs). Indeed 
for normal radiation, i. e. with 7 = 4, the amplitude of 
thermal density fluctuations on a given scale L can be 
expressed in terms of the entropy contained in a sphere 
of radius L (eq. (15.26) of [1]): 



(46) 



Thus, if entropy is produced at different rates on different 
scales the spectrum of fluctuations could be modified. In 
addition this issue is bound to interfere with any normal- 
ization condition for the amplitude. 



Harrisson-Zeldovich spectrum plus gravitational instabil- 
ity rather than a "proof" of the inflationary origin of the 
initial fluctuations. 
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VI. CONCLUSIONS 

We have studied the possibility of thermal fluctuations 
providing seeds for currently observed cosmic structures. 
Assuming thermal equilibrium, we have identified the 
necessary conditions under which these fluctuations are 
Gaussian and scale- invariant. We have shown that Gaus- 
sianity constraints are automatically satisfied in thermal 
equilibrium. The situation with the scale-invariance of 
the power spectrum is more problematic. In order to 
have a HZ spectrum in a universe that cools in time one 
needs significant modifications to the Stefan-Boltzmann 
law, namely, p oc T 7 with 7 < 1. We have shown that 
this condition is unlikely to be satisfied for radiation com- 
prised of Bose-Einstein particles. If we are prepared to 
assume that the radiation is not made of particles, then 
naturally a number of possibilities could open up. If over 
a certain period the Universe was warming up with time, 
then the condition for scale invariance becomes 7 > 1 and 
can, in principle, be achieved with a gas of Bose-Einstein 
particles with appropriate dispersion relations. 

An approach different from ours was taken in Ref. [30] . 
There, authors have considered thermal fluctuations in a 
system of particles with an attractive short range poten- 
tial and a repulsive 1/r 2 potential at large scales. They 
have shown that the resulting power spectrum of density 
fluctuations is scale-invariant on sufficiently large scales 
(small fcs). We find this direction of thought interesting 
and deserving further development. 

To conclude, we have discussed an alternative, "ther- 
mal" , way of obtaining scale-invariant initial spectrum 
of density fluctuations. Other ways to match the obser- 
vations without inflationary quantum fluctuations have 
been previously discussed, e. g. active seed models [31] 
with [32] or without VSL [33]. One may argue that 
some of these models are less "natural" than inflation. 
However, they make the point that the recent observa- 
tional victories in modern cosmology are a success of the 
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